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Abstract: The very notion of a current fluctuation is problematic in the quantum context. 
We study that problem in the context of nonequilibrium statistical mechanics, both in a 
microscopic setup and in a Markovian model. Our answer is based on a rigorous result that 
relates the weak coupling limit of fluctuations of reservoir observables under a global unitary 
evolution with the statistics of the so-called quantum trajectories. These quantum trajectories 
are frequently considered in the context of quantum optics, but they remain useful for more 
general nonequilibrium systems. In contrast with the approaches found in the literature, we 
do not assume that the system is continuously monitored. Instead, our starting point is a 
relatively realistic unitary dynamics of the full system 
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1 Introduction 

Certain aspects in the combination of nonequilibrium physics with quantum 
theory are often more problematic than their counterparts in nonequilibrium 
classical statistical mechanics. An important reason is that in statistical me- 
chanics one often starts from fluctuation theory and from estimates of statistical 
deviations. But in the quantum case, even for equilibrium statistical mechan- 
ics, there is no standard large deviation theory (although recently progress was 
made in [38l[3ll[22] and in the nonequilibrium case in [1]). For nonequilibrium 
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purposes one wants to go beyond first order perturbation theory around equi- 
librium, i.e. beyond covariance estimates. The question then emerges what one 
accepts as the definition of (also large) fluctuations of heat, work and currents. 

As it is often the case, such questions are more prone to confusion when 
working in the quantum formalism. Since recent developments in noncquilib- 
rium physics have focused on fluctuations (of entropy production), quantum 
analogues have been attempted by many different groups, etc. [3j>l G2E SS H21 

\m nam m m- 

Throughout this article, we choose the setup of a system connected to heat 
reservoirs and the fluctuations we study, are fluctuations of the heat currents. 

In Section [21 we present two possible approaches to fluctuations in a Hamil- 
tonian setup. These approaches have appeared repeatedly in the above-quoted 
articles. We remark that they are equivalent as far as the mean and the variance 
are concerned. 

In Section[3j we put the Hamiltonian description aside in favor of an effective 
model - a quantum master equation derived in the so-called weak coupling limit. 
In the framework of this master equation, we can again distinguish different 
approaches. One of these is the formalism of quantum trajectories, which is 
discussed in Section [H In Section El we combine the Hamiltonian description 
with the effective model. Our main result is contained in formula ()5. 1|) . It 
states that the fluctuations calculated by quantum trajectories are limits of the 
fluctuations in the Hamiltonian description. This result supplements the well- 
known derivation of the master equation, be it its rigorous form, as in [14] . or 
a more heuristic derivation, as in [T2] . 

One can further remark that both the fluctuations in the quantum trajec- 
tory picture [42] and the fluctuations in our Hamiltonian description satisfy the 
celebrated Gallavotti-Cohen fluctuation theorem, see 00] for a discussion. How- 
ever, stressing this point would be misleading because one does not need the 
weak-coupling limit to state the Gallavotti-Cohen fluctuation theorem. 

The work is mostly inspired by and based on [32 [TH] . Section E] presents a 
summary of mathematical arguments and states the main message of this paper 
as a theorem. 

2 What is a current fluctuation? 

2.1 Question 

Imagine several heat reservoirs Rfc, indexed by k € K . Each reservoir is in ther- 
mal equilibrium at inverse temperatures (3k and well separated from the others. 
All heat reservoirs are connected with a small system S through a coupling term 
proportional to a small coupling constant A. Formally, the composite system is 
described by a quantum Hamiltonian 

fceif fceif 

in which one assumes a clear separation between the reservoir Hamiltonians 
Hft k and the interaction i?s-R fc - The dynamics of the full system S + Rfc 
is given by the unitary group e~ ltHx on a Hilbert space of the form "system £g> 
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reservoirs." Imagine that the coupling between the system and the reservoirs 
starts at a certain initial moment. We take the initial state represented by a 
density matrix po of the form 



Po = Ps 



keK 



(2.2) 



where the states p k ,p k are equilibrium states at /3k on the fc'th reservoir, and ps 
is an arbitrary density matrix on S. 



We want to ask how much energy has flown out of/into the different reservoirs 
after some time t, and how this quantity fluctuates. In the classical setup, there 
is no ambiguity as to what that means: There one has a phase space X for the 
total system with a Hamiltonian flow ihi ( and the object of interest is the 
variable (function on X) 

x»H Rk (x t )-H Rk (x), (2.3) 

where Hn k now also represents the energy of the reservoir. Usually (but not 
necessarily), one starts from an initial distribution po, rather than from a fixed 
phase space point x S X, such that the above variables are actually random 
variables subject to some overall constraints like total energy conservation. A 
natural method to study the fluctuations of (|2.3p is to proceed via its charac- 
teristic function, for k € M} K \, with coefficients k^: 

J p a (dx) e-^eK ^{HK k (x t )-H Kk (x)) _ ( 2 4 ) 

Often, this formula is expressed in terms of a time- integrated current. 

We now switch back to the quantum case. It seems that there is more than 
one way to generalize (|2.4[) , as has been remarked by several authors. 



2.2 Answer 1 



The question amounts to choosing a quantization of the observable contained 
in (12. 4|) . Following the usual quantum dictionary one is tempted to introduce 
a "current operator" 

I k (t):=-iU* t [H x ,H Rk ]U t x , (2.5) 

where U x :— e~ ltHx is the dynamics generated by the total Hamiltonian H\ and 
-Hrj. is the free Hamiltonian of the fc'th reservoir only. Obviously, 



U x t H Rk U t x - H Kk = I dul k (u). 



(2.6) 



and one might set out to study fluctuations of the heat by considering fluctua- 
tions of the operator (|2.6p . This amounts to replacing (|2.4[) directly with 



Po 



(2.7) 



The expression (I2.7|) looks rather elegant but we do not know of any experiment 
or theoretical consideration where this quantity enters naturally. Observe for 
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example that U_ t H^ k £/ t does not in general commute with H^ k , hence their 
difference does not seem to be easily measurable. We therefore prefer a more 
operational definition that we present in the next section. Nevertheless, the 
current operator as defined through (|2.5p has its place in the literature, e.g. in 
the quantum formulation of the Green-Kubo relations j45[ [25] . As we remark 
in Section [2111 the definition (|2.5[) coincides in the first and second order with 
the definition which we present below. Hence, for the Green-Kubo relation, it 
does not matter which definition of current fluctuations one chooses. 



2.3 Answer 2 

A quantity that seems more relevant practically is the following: Assume for 
simplicity that (H^ k )k£K have discrete spectrum, indicating that we have not 
taken the thermodynamic limit and let x € X label a complete set of eigenvec- 
tors \x) of (H$i k )k£K with eigenvalues (H-R k )keK{x)- The corresponding spectral 
projections are denoted P x := |a;)(a;| and, by a slight abuse of notation, we use 
the same symbol P x to denote 1 ® P x , where 1 is the identity on Hs ■ Then we 
define the characteristic function as 

X (M,A,po) := ^[PyU^P xP0 P x U^ t P y ]e-' l ^^ K ^ H ^ v) - H ^). 

x,y£X 

(2.8) 

The idea behind this formula is clear: measure the reservoir energies (thereby 
projecting the reservoirs on the eigenstates x), at time s — switch on the 
interacting time evolution Ug, at time s = t switch the interaction off, and 
finally measure again the reservoir energies (projecting on the eigenstates y). 

We now use that the initial state po is diagonal in the basis |x) to rewrite 
([2~8|) as 

x(n,t,\,p Q ) = po \e~ i ^^^ H * k u t x e i ^>c^H IikU \] _ ^ 

Actually, we take the expression (|2.9p just as our starting point. Our main 
result is valid only after taking thermodynamic limit, in which the operators 
(i?R, fc )&gK have continuous spectrum and the definition of the state po has to be 
reconsidered. Nevertheless we will see further on that (|2.9I) has a well defined 
thermodynamic limit. 

Usually, when considering a system interacting with reservoirs, it is assumed 
that only system observables can be measured, since the reservoirs are very large 
and difficult to control. (|2.9j) does not follow this rule: it involves measuring 
reservoir observables Hn k . Note, however, that H^ k are reservoir observables of 
a special kind: they commute with and hence they are constants of motion 
for the free dynamics, which acts outside the time interval [0,t]. Measuring of 
i?R fc at times and t is conceivable even if the reservoirs are large (but finite), 
since to do this we have an infinite amount of time: s£]- oo, 0] for the initial 
measurement and s € [t,oo[ for the final one. Therefore, in our opinion, (|2.9|) 
can be viewed as measurable in realistic experiments, even though it involves 
reservoir observables. 

This approach to fluctuations was already used in [23] for fluctuations of 
heat, in [41] and very recently in [44] for fluctuations of work, and, most 
widespread, starting in [52" ] 155 ] . for fluctuations of charge transport. Note also 
the elegant approach to statistics of charge transport in [551 1Z] • 
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2.4 Comparison 



As we have seen, the question "what is a current fluctuation?" does not have 
a unique answer. In any case, both definitions of current fluctuations we have 
discussed above coincide in the first and second moment. In other words, the 
first and second derivatives with respect to k of expressions (|2.7p and (|2.9p 
coincide and are equal to 



Po 



I k (u)du 



Po 



I k (u)du / I k ,(u')du' 



(2.10) 



or, alternatively, to 

Po P-tH Rk U t -H Rh ], po [(U- t H Rk U t - H Rk ) {U. t H Rh , U t - H Rk , )] . (2.11) 

This can be easily checked (disregarding possible subtleties due to unbounded- 
ness of operators) by using that [H Rk , H Rk , ) = and that po is diagonal in H Rk , 
i.e. for all operators A: 

po[H Rk A]=p [AH Rh }. (2.12) 

Conclusion: If one is interested in second order fluctuations, for example 
Green-Kubo and Onsager relations, both generating functions (|2.7p and (|2.9p 
are equivalent. 

Note that the equality between (|2.7|) and (|2.9p up to second order depends 
crucially on the choice of the initial state. However, often in physics, one is 
interested in expressions which are independent of the initial state. 

Assume that the dynamics U-t ■ Ut admits a unique nonequilibrium steady 
state /9oc0- One expects that the covariance 



1 



lim -po 

t| + oo t 



(2.13) 



du{I k {u) - po [/*(«)]) / du'{I k ,{u')- po[hi{u')]) 
'o Jo 

which is a combination (see also (|2.15p ) of the two expressions in (|2.10p for t j oo, 
is independent of the initial state and in particular equal to the correlation 
function 



du p c 



[h{u) - Poo [4(o)] ) O(o) - Poo [i y (o)] ; 



(2.14) 



Using the above reasoning and standard manipulations of the cumulant gener- 
ating function, one easily checks that (|2.14p equals 



lim 



if- 



if+oo dK k dn k i t 



7 log 



Expression 



(2.15) 



When poo is an equilibrium state, hence all temperatures equal, then expres- 
sion (12.14p with poo [4(0)] = features in the Green-Kubo relation, as stated 
rigorously in [25] and specifically for quasi- free systems in [3] . 

Having chosen a definition of a current fluctuation in a microscopic (Hamilto- 
nian) description, we set out to consider an effective model, arising by a certain 
consistent approximation. The model we will be dealing with in the present 
paper is that of the quantum master equation which often arises in the so-called 
weak coupling limit. We give some reminders in the next section, and we con- 
tinue the answer to our question in Section [3J In Appendix B, we outline the 
weak coupling limit of both generating functions (|2.7p and (|2.9p . 

4 Obviously, one needs the thermodynamical limit for this assumption to be realistic 
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3 Weak coupling limit 



One of the aims of the present paper is to supplement the user's manual for one 
of the best known effective equations in the physics of open quantum systems: 
the master equation for the long-time evolution of the density matrix of a small 
system with discrete spectrum in contact with reservoirs. It is widely accepted 
that master equations gives a good description of the degrees of freedom of a 
small system in certain limiting regimes. In the physics literature, this limiting 
regime is characterized by the Born-Markov approximation and the rotating 
wave approximation, see for a review. Another name, common especially 

in the mathematical physics literature, is the weak coupling limit, which makes 
these two approximations exact. It goes back to |23j and was made precise in 
[14] . An interesting review is contained in [3D]. We start with a formal sketch 
of the usual set-up. 

The small system is modeled by a finite-dimensional Hilbert space S and a 
Hamiltonian H§ — a certain Hermitian matrix. It interacts with an environment, 
possibly containing several reservoirs indexed by k £ K. Let us think of each 
reservoir as an assembly of free oscillators. We use the well-known notation 

dT(h k ) = dqh k (q)a* k (q)a k (q), (3.1) 

where h k is a dispersion function on R d , the so called one-particle Hamiltonian 
on the one-particle Hilbert space t) k := L 2 (R d ) and dT(h k ) acts on r s (f)fc), the 
bosonic Fock space corresponding to \) k . 

The coupling between the system and the environment is linear in the sense 
that 

H S - Rk =[ dqV k ^[f k (q)a* k (q)+JMa k (q)], (3.2) 

where f k G \) k and V k are self-adjoint operators on Tig. The total Hamiltonian 
is formally given by 

H X =H S ®1+ ^l®dr(/i fe ) + A^ffs-R fe . (3.3) 
fee/f fceif 

on the Hilbert space Tis <B> [® k £KF s (t) k )]. Observe that the prefactor A measures 
the interaction strength. 

The Hamiltonian H\ generates a quantum evolution = e~ lrHx . The 
weak coupling limit concerns the convergence of the reduced dynamics on the 
small system. The coupling A J, gets very weak as the time r = t/X 2 goes to 
infinity. A well-known theorem by Davies [T3] states the following result (which 
is written here in a formal way, one actually needs the framework of Section [H] 
or a limiting procedure like in Appendix A to give it a precise meaning): 

lim a, [Ul H U\-H(S®l)Ul- H U\- H ] = p s [e tc S] (3.4) 

for matrices S on Tls- The initial state po = ps <£> [® k ^KP k ,i3 k ] is the product 
of an arbitrary state (density matrix) ps on Tig and of thermal states p k ./3 k 
at inverse temperatures k in the respective reservoirs H k . The superoperator 
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(acting on matrices) C is the generator of a completely positive dynamics ob- 
tained in the weak coupling limit. It can be written in the so called Lindblad 
form [S3], see (|4.ip below. There are of course technical conditions for (|3 .4[) 
to be true, which we skip here. One should realize that (13. 4|l is a non-trivial 
statement. It specifies conditions under which the (reduced) evolution on the 
small system S gets autonomous. The result is however somewhat rough, at 
least for nonequilibrium practice, as the resulting process is just a jump process 
between energy levels of the system. In particular, it does not allow us to track 
the interactions with a given reservoir. What follows is a way to remedy that, 
at least on a phenomenological level. 

A first observation (which can be checked from the explicit construction 
given in Sections [H and [5]) is that the generator splits naturally as 

k&K 

where each Ck can be defined as the object that would emerge in the weak 
coupling limit from the microscopic Hamiltonian by cutting the interaction with 
all spaces ()& except for k' = k. 

Let us see how the decomposition (|3.5p can inspire us further to answer 
the question of Section 2 in the case of master equations. One can define the 
time-evolved current operators 

J fc (i) =e tc (£ k (H s )), (3.6) 

(much in the spirit of (|2.5[) ). where we recall that H§ is the Hamiltonian of the 
small system. 

Now, one might conjecture that the characteristic function 



Pa 



-i£* K k Jo duj k (u) 



(3.7) 



is the limit of the correlation function (|2.7p in the weak coupling limit. That is 
actually not correct! Neither is it correct that (|3.7[) is the limit of (|2.9[) . (|3.7[) . 
however, coincides with the limits of both (|2.7p and (|2.9p as far as the mean 
current (first moment, i.e. first derivative in k) is concerned. 

Nevertheless, in [30], one starts from the current operators (|3.6[l and one 
obtains the correct Green-Kubo relations. This is due to the fact that one does 
not calculate the variance of the current via the characteristic function (|3.7p , 
but instead, one starts from the current-current correlation function which is an 
analogue of (|2.14|) . (See [32] for a more general treatment of the Green-Kubo 
relations in the weak coupling limit) 

In the next section, we present another (better) way to identify the fluctua- 
tions in the weak coupling limit. 



4 Quantum trajectories 



Let us look a bit closer at each of the components C k of the weak coupling 
generator. They are given by 
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where c(w, k) are positive constants, Ek are effective Hamiltonians, sometimes 
cailed Lamb-shifts, and 

V u , k := J2 PeVkPe', (4.2) 

e,e' £spHs ,e — e'—uj 

where now P e are spectrai projections of Hs corresponding to the eigenvalue e. 
The summation in (|4.ip is over all differences of eigenvalues of Hs (or equiva- 
lently, over all eigenvalues of [Hs, •]) and in what follows we write simply ^ w 

for E w esp([ff 3) .]) and Efc for Efcex- 0ne now decomposes £ k = £° k + £ w Ju,,k 
with 

£2(5) = i[S fc , 5] - ~ c ( w > *){K,,fcV* fc , 5}, X, fe (5) = c(w, fc)V w ,fcSV* fc , 

(4.3) 

where j7L,fc is called a jump operator. It is important to keep in mind that such 
a splitting C k = C k +Juj,k is not uniquely given by the generator £/., instead we 
have used information about the operators Hs and V k to define this splitting. 
(See j!2j for extensive comments on this non- uniqueness). 
The final unraveling is written as 

C = C + ^J u , k , (4.4) 

uj,k 

where Co = J^k^k- We now introduce completely positive operations Wt{d) 
and Aj given by 

W t (a) := klJ UlM Kl- tl ■ ■■^l-t n . 1 J^,kA° t -t n , A° = e t£ °, (4.5) 

for a "trajectory" a that labels all the jump times and actions in Wt{u)'- 

a= (ti,kx,u 1 ;t2,k2,0J2;---;t n ,k n ,Un), < tx < .. . < t n < t. (4.6) 

Via the (norm convergent) Dyson expansion corresponding to the splitting (14. 4|) . 
we have 

e tc = J daWt(a), (4.7) 
where the integral over a is the abbreviation of the following expression: 

v ft pt n f 

12 12 12 / dtn / dtn - 1 ■■ dti = : / da - 

n=0k 1 ,...,k n u; 1 ,...,u n J ° J ° J J 

The idea is now to interpret each contribution to the sum and integral as a 
quantum trajectory. In that philosophy, the map 

S^W t (o-){S) (4.8) 

gives the (unnormalized) evolution of the system, conditioned on the trajectory 
a. This conditioning usually means that certain measurement outcomes were 
obtained and that these outcomes are represented by cr. 

Taking this idea just one step further, one can obtain statistics of mea- 
surement outcomes (or, in our case, currents). Define the following probability 
distribution on all possible a: 

dP* s (<7):=Tr[p s W i (cr)(l)]dcr, (4.9) 
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and introduce the energy counting numbers n\. 



nl{cr) :=^^,ftWl( fi <t)( ff KW» ( 4 - 10 ) 

i 

where the index i runs over all jumps present in a (i.e. from 1 to n in (|4.6p ) and 
l(ti<t) is the indicator function of the event that in a, the i'th jump occurs before 
time t. The distribution of the random variables n\, is inherited from (|4.9p and 
will be used throughout. The characteristic function of the joint distribution on 



L k£K 

Xw.c.(M,Ps) := / dP*» [e-^^^tW] . (4.11) 

Note that (|4.11[) characterizes the full distribution of n\ eK . The study of fluc- 
tuations indeed amounts to more than characterizing covariances, as useful as 
that may be within linear response theory. 

A further point concerns the "classical" nature of the variables n\ (a) and the 
use of standard probability theory. But that is exactly the point of the present 
paper: these variables characterize the "quantum" fluctuations, see further in 

GB. 

In the quantum optics literature, this procedure of "unraveling master equa- 
tions into trajectories" is generally accepted, see e.g. [13], [12] and the recent 
review [TIT] , whereas the first source of unraveling is probably in [13] . 

From a more fundamental point of view, one could ask how quantum trajec- 
tories and, more specifically, the variables n\ emerge from microscopic dynamics. 
The usual justification of quantum unravelings found in the literature supposes 
that the system is described by a so-called quantum Langevin dynamics (the 
solution of a quantum stochastic differential equation (QSDE)). Even though it 
is a unitary dynamics, it can be viewed only as a very approximate description 
of realistic quantum systems quite far from "first principles" . To justify the use 
of QSDE it is usually assumed that the unitary evolution is interrupted by mea- 
surements, which, in the limit of very short times between measurements, yields 
a quantum Langevin equation (originally introduced by [24j in a mathematical 
framework, see [3] for recent developments and |21] for a physical point of view.) 
In many cases it amounts to supplementing the quantum evolution of the small 
system with a stochastic evolution of classical variables in the environment. 

The quantum stochastic differential equation (QSDE) or stochastic Schrodinger 
evolutions, and their solutions, have been studied by many authors. The ran- 
dom variables ri^ correspond there to the fluctuations of (linear combinations 
of) number operators in the reservoir spaces of the QSDE. (see e.g. [H] or [tJJ) 
We will not elaborate on this point, since it is not our central subject (See 
however Appendix B). 

Lately, a new class of models was considered, which give also rise to QSDE's. 
These are the so-called "repeated interaction models" , where, instead of mea- 
suring the reservoir continuously, one refreshes it, see [5] [5]. 
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5 Connecting unravelings with the Hamiltonian 
dynamics 

Our paper describes an alternative justification of quantum unravcllings that we 
believe is in many situations more satisfactory and does not explicitly introduce 
continual measurements, stochasticity or "refreshing of reservoirs" . We start 
from a class of dynamics described in Section [3j which are often viewed as a 
relatively adequate description of realistic quantum systems. We prove that, 
after first applying the thermodynamic limit and then taking the weak coupling 
limit, the quantities (|2.9p converge to quantum unravellings -expressions of the 
form (|4.1ip . This idea, to our knowledge, is present in the literature only in a 
heuristic form. What we describe is a rigorous result. 

Our proof of this result involves two steps. The first step is the theorem 
about the extended weak coupling limit obtained by two of us [18) . which says 
that the microscopic dynamics converges in a certain sense to an appropriate 
quantum Langevin dynamics. The second step (which is well known in the 
literature) goes from the quantum Langevin dynamics to quantum unravellings. 

There exists actually an alternative proof of our result that goes directly from 
the microscopic dynamics to quantum unravellings, without passing through a 
quantum Langevin dynamics. We will indicate it briefly in Remark [ 

Recall that 



H Rk = dT(h k ) 



dqh k (q)a* k (q)a k (q) 



is the (second quantized) k— th reservoir Hamiltonian, cfr. (13. Remember 
also the characteristic function Xw.c.(^, t, ps) from (14. lip . Here comes the main 
result of the paper. 

Under standard conditions of the weak coupling limit as in (|3.4p . and with 
the same remark about precision as in (|3.4|) . 



Xw.c.(k, t,p s ) = lim />o 



(5.1) 



The right-hand side is of course a general instance of the weak coupling limit 
of (|2.9p . This proves that quantum trajectories provide nonequilibrium fluctu- 
ations of the time-integrated heat dissipated in a given reservoir. In particular, 
the properties of the distribution of n k , such as these related to fluctuation 
theories or to Green-Kubo relations, are related to the quantum fluctuations 
of energy currents in the sense of (|5.ip and the arguments of Section [5J An 
extensive study of the distribution associated to ()3.7p was performed in [32] ■ 



6 Mathematical statement of the results 

The main result of the paper, formula (|5.ip . is a consequence of the more general 
and abstract results proven in [18] . For the convenience of the reader, we list 
some simple assumptions which allow to establish f|5. 1|) and we specify what is 
the exact definition of the quantities appearing on the RHS (|5.ip . For simplicity 
of presentation, we choose a special form for the one-particle Hamiltonians h k 
and we restrict the physical dimension of the reservoir space d to d = 1. 
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6.1 Formal Hamiltonian 

Recall the formal Hamiltonian from (|3.3[) : 

Hx = H s +J2 H ^ k + AV* ® (a* k (f k ) + a k (f k )), (6.1) 
fee k 

where now 

• H s = H£,V k = V k * e B(Hs) with dining < oo; 

• Hn k = dT(h k ), where h k are the one-particle Hamiltonians on the Hilbert 
spaces f)fc = L 2 (R + ) acting as 

(h k g)(x) = xg(x); (6.2) 

• fk S f)fc are coupling functions; 

• a- k (f k )/ci k (f k ) are creation/annihilation operators on the bosonic Fock 
space r s (f) fe ); 

6.2 Effective master equation 

Given the information of Section 16.11 and the inverse temperatures (3 k , k S 
K, one can construct the weak-coupling generator £, which was introduced in 
Section [3] with unspecified parameters c(uj,k) and E k . Define the functions 




(6.3) 



The exact expressions for the parameters c(w, k) and E k are (see e.g. |18j ) 

c(^,fc) = ^IJfMI 2 , (6-4) 

Z7T 

Sfe = E F "V^/ die_iU ' f dxe- itx \f^(x)\ 2 . (6.5) 
6.3 Coupling to thermal reservoirs 

One of the subtle points of quantum statistical physics is how to describe in- 
finitely extended bosonic reservoirs at positive temperatures. Strictly speaking, 
the Hamiltonian H\ defined in (|6.ip describes the reservoirs only at zero tem- 
perature - but we are interested in the case of an arbitrary temperature. (This 
is why H\ was called the "formal Hamiltonian"). We need the state po in ex- 
pressions p.4[) and (|5.ip to be a thermal state (which cannot be represented by 
a density matrix in infinite volume). In expression (|3.41 15. l|) . we pretended that 
this state can be defined on B(li.s® [®kEK~T s 0)k)]), or at least on its sufficiently 
large subalgebra preserved by the dynamics, but this could be problematic. 
(This would however be a good approach for fermions!) 

Fortunately, there exists a formalism that allows us to describe a system 
interacting with reservoirs at positive temperatures in the thermodynamic limit 
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rigorously. This formalism was used in standard works like [57J HS1 IB]- The 
relevance of this construction has been argued e.g. in [T^l [TTJ- Here we just 
present how one should modify the dynamics of the coupled system and we 
point to Appendix A for a justification. 

One of the ingredients of this formalism are the so-called Araki- Woods rep- 
resentations of the CCR. To introduce them one needs to enlarge the Hilbert 
space. The enlarged Hilbert space is 

H :=H s ®T s (® keK (t) k ® t) fe )). (6.6) 

We define the free Liouvillian of the fc'th reservoir on T s (t) k © f)fc) as 

L Rk =6T(h k ®{-h k )). (6.7) 

From now on, it will be convenient to identify t) k ffi with L 2 (R) such that 
the one-particle operator h k ffi (—h k ) acts by multiplication with x € R. The 
generator of the dynamics is chosen as the so-called semi- Liouvillian (see e.g. 
[16j for explanations about the terminology) and it equals 

La = H S + X Vk ® + MIL 3 ")) + J2 L ^- ( 6 - 8 ) 

keK keK 

This is a formal expression, but one can easily construct the operator L\ rigor- 
ously (see [IB])- Finally, let f2 stand for the vacuum vector in T s ((BkeK(^)k®i)k)) 
and define the vacuum state 



Vac[-] = (ft, -ft). (6.9) 
We can now define our object of interest: 

X[k, t, A, ps) = (p s «)Vac) L" i2: ^« Kfci ^e- itiA e iE '=^ KfcLR ^e itiA j . (6.10) 
Remark that (pHX)]) arises from (|2.9p . that is 



by replacing 



® PKfJ k 

keK 



by Vac, H^ k by L^ k and H\ by L\. One can check 



that in finite volume both expressions coincide - in particular, the positive 
temperatures of the reservoirs have been incorporated directly into the functions 
fP k . In fact, the vacuum state Vac represents the product of thermal states on 
the appropriate algebra of observables. More details and explanations can be 
found e.g. in [T5] , 

In Appendix A, a limiting procedure which constructs expression (|6.10p via 
finite-volume approximations is explained. 



6.4 Main result 

We are ready to state rigorously our main result 

Theorem 6.1. Let Xw.c. ( K : t, ps) be defined via (|4.1ip as in Section^and with 
parameters \6.4\6- 5\) . Assume that 
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1) The reservoir correlation functions are integrable 



At 



dx\f£ k (x)\ 2 e- itx <oo. (6.12) 



2) The functions f^ k are continuous in a neighbourhood of (spHg — sp-ffs). 
Then for all k G R'"^ , 

limx('«, A~ 2 t, A, ps) = Xw.c.(K,t,p s ). (6.13) 

AJ.0 

This theorem is a simple consequence of Theorem 5.7 in [15] . Essentially, one 
proves that e _l4A La converges to a quantum Langevin dynamics and operators 
like e _IKfcA converge to e~ 1KkNk , where Nk is an appropriate linear combi- 

nation of number operators (see also Appendix B). The link with unravclings 
of master equations belongs then to standard knowledge in quantum stochastic 
calculus (e.g. to be found in various forms in [4"2l [TT1 15]). 

Remark 6.2. The assumptions contain a mild infrared regularity requirement, 
since the assumption about f^ k implies that x i— > x^ 1 ! 2 fk(x) G L 2 (R + ). 

Remark 6.3. Note that one would like to strengthen the statement of Theorem 
1 6. II For example, one would like to have convergence of the derivatives in k. This 
is possible under stronger regularity assumptions on fj~. One can also prove a 
version of Theorem \6.1\ allowing for complex Kk, see [39]. 



Remark 6.4. Although the extended weak coupling limit in I7<SJ / gives valuable 
insight into the limit of fluctuations (see Appendix B), one does not really need 
it to prove Theorem \6.1\ One can also rewrite %(k, t, A, po) a s 

X(k, t, A, po) = (ps ® Vac) [ e - itL ^/V tL ^-^ 2 ] . (6.14) 

where 

Lx, K = H s +Y,Ln k (6.15) 

keK 

+A V k ® (<(e- iKfc( ^ e ~ hfc) /f fe ) + a fc (e- Efc( ^ e ~^ ) /f' i )) , 

which, at least for 3k = 0, reduces technically to the usual derivation of the 
master equation. 

APPENDIX A 

We justify (|6.10p from the physical point of view. We do that by arguing that 
(|6.10p is the thermodynamic limit of finite volume versions of (|6.1ip . 

Let for each n G N, f)k,n be finite-dimensional Hilbert spaces with one- 
particle Hamiltonians hk, n and coupling functions fk, n G f)k,n- Define the finite 
volume evolution reservoir Hamiltonians Hn htn and full Hamiltonian H\ n by 

#Rk,n = dr(/ifc jW ) 

H x ,n = Hs+Yj HR *>" + X Ys Vk ®( a *(f k ^ +a (f k > n ^ 

k£K keK 
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and the finite volume thermal states on £>(r s (f)fc. n )) as 

Pk A , n [■} = Tr [e" &ffR *- •] , Zk,n(Pk) = Tr [e-A>***.»] (A-l) 

Zk,n{Pk) 

Assume that 

1) There is a C < oo such that for alii G R, A; G K and n G N, 

\pk A ,n [a*(e- it ^."/ fc ^)a(/ fc ,„)] | < C 

2) For each teR and k e K, 

lim p fcA ,„ [a*(e- it ^.»/ fc , r[ )a(/ fc ,„)] = (A, (e^* h * - l) _1 e- it ' , */fc>6» 

nf oo 

lim p fc)A)n [a(/ fc ,„)a*(e- it ^-/ fc ,„)] = (A, (1 - e"^^)-^"^ A)„ fc 

n]'oo 

The notation on the RHS was introduced in Section [3] 
Then the expressions 

X„(M,A,Ps) := (ps®[ ® Pfc,ft,n]) f e - i S*' , * HR ».»e- i " r *'»e i S* ,e * HR »-»e i " r *'" 

(A-2) 

converge as ii f oo for all i and A. This is checked by writing a Dyson ex- 
pansion for (|A-2[) . treating the terms in A as a perturbation. The dominated 
convergence theorem can be applied since term-by-term convergence is implied 
by Assumption (2) above and a dominating bound follows from Assumption (1) 
above. 

The connection with the setup in Section El is that, under the above assump- 
tions 

lim Xu(k, t, A, p s ) = t, A, p s ) (A-3) 

njoo 



APPENDIX B 



In this appendix, we look more closely at what the "extended weak coupling 
limit" (as presented extensively in [18]) can tell us about the different fluctuation 
formulas that were proposed in Section [2] 

Assume the notation introduced in Sections 2 and 3. The result in 1181 states 



that the unitary dynamics U^ 2t on Hs <8>fc r s (f)fc) converges (in an appropriate 
sense) to a new unitary dynamics Ut on Ws<S)fc,a>r s (f)fc ia ,), where f)k tU1 are modified 
one-particle spaces. The dynamics Ut, which is the solution of a Quantum 
Stochastic Differential Equation, is extensively discussed in [4"2l I18j . 

An important observation is the emergence of effective reservoir energy op- 
erators; 



N k := 



(B-l) 



where l W) fc is the projector on i)k,u- One sees hence that the effective reservoir 
energy operator is more like a number operator. Its plays a role in the limits of 
respectively (|2.7[) and (|2.9j) . One has, of course again under technical conditions, 



Po 



A\0 



PO 



-iE feeK K k (U- t N k U- t -N k ) 



(B-2) 
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and 



(B-3) 

where po is a state which coincides with po on Tig and which represents the 
thermal reservoir states on (£>kF s (t)k)- Although it is not obvious from these 
expressions, the expression (|B-3|) coincides with (|4.11[) . this is the well-known 
connection between unravelings and quantum stochastic differential equations. 
To check that (|B-2[) and (|B-3|) coincide up to second order in k, it suffices to 
know that 

Pa[AN k ] = p [N k A] =0 (B-4) 
for all A, operators on Hs ®k,u r s (f)fc iW ). 
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